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We consider a lattice-inspired random matrix model for the QCD chiral phase transition at finite 
chemical potential. Useful features of the usual RMM for QCD at finite ^ are reobtained, some 
being brought closer to their lattice equivalent. The simple physical requirement of a vanishing 
quark number density in the broken phase is fulfilled in the limit of a large number of timeslices. 

It is argued that the suppression of the partition function at nonzero /i in the broken phase, seen in 
the usual RMM, is possibly present in lattice simulations and is simply a result of the discretization 
-_- ■ in time. 
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^ ^ The random matrix model (RMM) for chiral symmetry breaking in QCD at finite chemical potential has been 

^ introduced by Stephanov Significant qualitative insight has resulted from employing the RMM as a schematic 
model of [lattice] QCD at finite chemical potential The failure of the quenched approximation in lattice QCD at 
' finite /i was explained in this context by the non-analyticity of the Nf — > limit |l) . An extension of the RMM to 
^ . include both temperature and chemical potential offers a simple description of the tricritical point in QCD Q . 

It was noted early on that the partition function of the RMM at finite n is suppressed in the low density phase [|j , 

■ specifically, it depends on /i like ^ where N is the size of the random matrix. The suppression was assigned to 
averaging over the complex phase of the fermion determinant. The study of the Glasgow method Q applied to the 
same model Q showed that the method requires exponentially large ensembles for convergence. This property was 

■ directly traced to the smallness of the partition function close to the critical chemical potential. 

\ However, the decrease of the partition function with fi below the critical chemical potential in unphysical. One 

would expect that in the low density phase there is no /i dependence so that the number density (n) = df^ In Z is 

identically zero. In this paper we propose a different way to introduce the chemical potential dependence in a RMM, 

JL i more similar to lattice calculations. 
S-H, 

. In the next section we introduce the model and work out the effective theory. The model copies the time-structure of 
^ the lattice Dirac operator at finite chemical potential. Using standard RMT methods and some algebraic manipulations 
. ^ . we reduce the initial integral over an N x N matrix to one over an Nf x Nf matrix, where Nf is the number of quark 
flavors. Sec. Ill is devoted to the analysis of the partition function for one flavor, using the saddle point approximation. 
H ] We identify the phase structure and show the central result of this paper, the /^-independence of the low density phase 
. . 1 in the limit of a large number of timeslices. In Sec. IV we show numerical results on the eigenvalue distribution of the 
model Dirac operator and the zeros of the partition function. The eigenvalues show a strong qualitative similarity 
to actual lattice results. The zeros of the partition function trace out the boundary between the two phases in the 
complex plane of the parameters, reinforcing our saddle point analysis. We summarize and discuss our findings in the 
final section. The Appendix contains technical details used in the derivation of the effective partition function. 

II. MODEL 
A. Definition 

Consider a partition function of the generic form 

-wtrcct f m® In igC -f W{^J) ^ 



Z.K,) = /PCe— det( ,^ct:V(-,) 'm^ir ) ' 
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which is dictated by the chiral structure, mass (m) and chemical potential {fi) dependence of the QCD partition 
function. The integration is over an N x N complex matrix C, where VC is the Haar measure. The size parameter 
g can be seen as a measure of the strength of the interaction. 

The usual choice for the /i-dependence, W{^) — fj, is also the simplest one [|l|,||J^. We wish to build a model whose 
chemical potential dependence mimics more closely that of the lattice, where the chemical potential enters through 
factors of e*^ multiplying the forward and backward links 0. Therefore we imagine that we have Nt timeslices with 
Ns points in each of them, so that N = NtNs- In the absence of a temperature dependence, we choose 

Win) = (e'^SN, ® In^ - e-^S?;^ ® 1^,) , (2) 

where Sn is the forward step matrix of size N. This way each block of size Ns is coupled to the next and to the 
previous block through a factor e^^, respectively. 



B. Reducing the partition function 



We start by performing a standard Hubbard-Stratonowitz transformation. This leads to replacing the integration 
over the complex N x N matrices C with integration over an Nf x Nf complex matrix a. The size of the matrix 

under the determinant can be reduced by 2 via the identity det ^ ^ ^ ^ = dct{AD — BD^^CD) and taking into 

account that the resulting upper and lower right blocks commute, we obtain 



Vae' 



-Nticrcr^ 



det 



(M + ga) {M + gaf 1 at, - 1 at^ «) (e^SjVt - e ^S^^ 



(3) 



Note that the two remaining terms under the determinant may be diagonalized simultaneously. In the Nf > 1 case, 
one may calculate the exact partition function by switching to the distribution of the eigenvalues of a. In the following 
we will simply assume equal masses and a diagonal, effectively leaving us with one flavor. We also assume that Nt is 



even, Nt — '2.N[, so the square of the step matrix breaks up into 2x2 blocks, 
reduce the size of the matrices by two. Our partition function is now 



£aae~^""' det | [(to + ga){m + ga*) + 2]® Iat/ 



= StVj' ® I2, and so we can again 

(4) 



The matrix under the determinant is of the form considered in the Appendix, 



QL{a,x) 



I a 



X 

a ) 



(5) 



with L = Nf, a = (m + g(j){m + gcr*) + 2, and x — — e^^. The determinant is then 

dct(g) = Ai(a)^ + A2(a)^-(-)^(a;^+x-^) . 



(6) 



where Ai2(a) are the two solutions of A'^ — aA + 1 = 0. For simplicity we will assume that N't is even. By convention, 
if they are real, we deflne A = Ai to be the one larger than 1 in absolute value. With all this, we have 



fc=i 



-N 



2N,Nf 



(7) 



where N = NtNs, Nt = 27V/, a ^ {m + ga){m + ga) + 2. 
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III. PHASE STRUCTURE 



Let us consider the case of one flavor, Nf = 1. The partition function depends on two thermodynamic parameters, 
the quark mass m and the quark chemical potential, /i. The derivatives with respect to these two quantities are pro- 
portional to the chiral condensate x — -^9 In (w, /z)/9to and the quark number density n = j^dln Zff{m, fi)/d^. 
By changing variables to r = cr + m it is easy to see that x — (25R(cr)). For positive m, the integral over the phase of 
T has a saddle point at arg(T) = 0, so the integral over (Pt turns into a radial one. With a further change of variables 
to u = T^, we obtain 



du ■ 



(8) 



Let us first focus on the to ^ behavior, by setting m = in (H) and bearing in mind that % = 2{^pu). We may 
write the partition function as 

oo 

Z(/Lt) - y"(iwe-^" (A(a)^'' +A(a)-^'' -2cosh(2/LtAr,'))^^' , (9) 



where now a — + 2 and N — 2NsNl.. What is new about our partition function in comparison with the usual RMM 
is that here we have two potentially large numbers, N^. and N!,. Since Ns is a large number (in lattice simulations it 
is of the order of N^, e.g., 500 for Nt =8), it is justified to approximate Z by the saddle point result corresponding 
to the absolute maximum of the function $, 

A(ug2 + 2)^'' + X{ug^ + 2)-^'' - 2 cosh{2fiN^)) j''^^ = g, N[, ^if'^' . (10) 

In actual lattice simulations the number of timeslices Nt = 2Nl, is not very large, for instance, in A^t = 2Nf — 8. 

For small Nt values out model can be worked out in the traditional manner. The integrand for N^ = 2 and g = 1 
is rewritten as 

{e-^"((u + 2f -4 cosh(2/i)2)}''^^ =$(«,<;, 2, . (11) 

The stationary points of the function are at ui2 = (^—3 ± + 8cosh(4/i)^ /2. Since the integration is from u — 0, the 

negative stationary point does not matter, but one needs to take into consideration the endpoint, u — 0. Altogether, 
we find two phases, one with x = {^/u) = for large /i and one with x for small fi. They are separated by a first 
order phase transition at /ic — 0.1966, when <i>(0, g = I, iV/ = 2, n^) = ^{ui,g = 1, Nt = 2, fic)- Similarly to the usual 
random matrix model at finite /i, the quark density starts off from zero ai /j, — 0, decreases until fic, where it has a 
finite upward jump, and then continues to decrease. 

For Nl larger than 4 a closed form solution is not possible. However, the general features of the function $ for 
M > are similar, in the sense that <&(m, Nj., fi = 0) starts off at zero, increases and has exactly one maximum, after 
which it asymptotically decreases again to zero. For > 0, the value at u = is negative, but the initial increase and 
the unique maximum are preserved for moderate values. 

We plot u, the location u of the maximum of $(u), as a function of /i for various N, in Figure |l|. We can see how 
this quantity, which is directly related to the chiral condensate in the broken phase, x — 2Vu, becomes increasingly 
independent of /i as Nt increases. In figure |^ we plot the critical chemical potential /ic for various values of Nt — 2N[. 
Again, fidNf ) quickly converge towards a limiting value. 
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FIG. 1. The location u of the maximum of 4>(it,g = with respect to u as a function of /i for various numbers of 

timeslices, Nt = 27V(. As Nt increases, the /i-dependence becomes negligible, leading to a low density phase with identically 
zero quark number. 
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FIG. 2. The value of /.ic for various numbers of timesteps Nt ~ 2Nl, at g — 1, m — 0. We define Hc as the value of mu at 
which the magnitude of <1> at u = equals the value of "I> at its maximum. 

Assuming Nt large, things become very clear. The quantities \{ug^ + 2)^^, solutions of a quadratic equation, are 
always real (for positive u), and one of them, A = Ai is larger than 1, and increases with u. Therefore the function 
e~^t"A(ug^ + 2)^t has a unique maximum (for positive u) at u = \Jl + A/g^ — {u — 0.236... for g — 1), 

corresponding to A = 5^/2 + ^71 + (?'*/4 (A = 1.618 . . . for g = 1). For a significantly large N'^ this maximum is not 
disturbed by the 1/A term. In order for the cosh(2/iAf/) term to matter, exp(2/i) has to be comparable to A, which 
corresponds to fj, around /ii = 0.2406 . . .. However, this does not come into play since the maximum at u competes 
with the absolute value of $ at u = 0, and the approximate condition for the equality of the $ values is 

exp{~2fic) — exp{~u)\ , (12) 

which gives fic = 0.1226 . . . < /zi . Therefore the maximum at u ceases to dominate before its location is significantly 
influenced by u. The transition between the maximum at u and the value at u = becomes sharper as Nt increases. 

In Figures and ^ we plot the saddle point result for the partition function and the two main observables, the 
number density and the chiral condensate, for a large number of spatial points [Ns » 1), and for various N't values. 
In Figure ^ we plot the logarithm of the absolute value of the partition function, normalized to its value at /i = 0, 
and divided by = NgNt = 2NsNt. As we have discussed, the /i-dependence for small Nt values is very similar to 
that seen in the usual random matrix model. The partition function decreases in magnitude with increasing fi up to 
the critical chemical potential. As a matter of fact, this suppression persists for all finite Nt but it drops quickly in 
magnitude as Nt increases, so that for Nt = 40 the partition function is practically flat between and ^c- The same 
trend can be seen in the plot of the number density ^. The decrease of the partition function is translated into a 
negative number density, which is certainly unphysical. This feature was present in the usual random matrix models, 
and was discarded as an artifact of the model. Finally, the chiral condensate is also sensitive to even in the low 
density (or broken) phase. This corresponds directly to the fact that the maximum of $ moves as /x varies. 
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FIG. 3. The absolute value of the logarithm of partition function at m = 0, normalized to its value 2% — Zn^P^ ~ 0) and 
divided by TV = NsNt — 2NsNl, in the saddle point approximation with respect to Ns, for various values of Nl = Nt/2. The 
value Nt — 2N't — 8 corresponds to usual lattice simulations. Notice that for larger Nt values the partition function becomes 
practically flat below fic- 
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FIG. 4. The value of the number density versus the chemical potential, at m = 0, corresponding to the partition function 
plotted in Fig.|[ 
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FIG. 5. The value of the chiral condensate versus the chemical potential, at m = 0, corresponding to the partition function 
plotted in Fig.^ The number density and the chiral condensate mutually exclude each other. 

An important conclusion we may draw from our analysis is that the physical requirement of /.t-independence of 
the low density phase is ensured by a large Nt, i.e., a fine enough time discretization. Another interesting point is 
that the suppression is always present to some extent. This is probably so for Glasgow simulations at finite /i, where 
the number of timeslices used was as low as Nt = 8. But even for higher numbers of timeslices, the suppression is 
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enhanced by a factor of NtNg. Although the suppression decreases significantly faster than l/iVj, the factor A''s is 
typically very large in lattice simulations, since it is of the order N^. Therefore it is not unlikely that the unphysical 
artifact seen in the usual RMM, leading to a negative number density, is also present in traditional lattice simulations. 
In our analysis of the Glasgow method, we have found a connection between the failure of the Glasgow method and 
the suppression of the partition function. 

For m 7^ 0, g 1 the analysis is similar. The iVt — > oo limit leads to a sharper transition and increasingly fi- 
independent broken phase. The mass dependence of the partition function is similar to that in the old model. For a 
positive mass, the exponential in $ has a maximum at uq = /g'^, which competes now with the maximum at u. 
The location of u is also shifted upwards and so is the size of the maximum, while the magnitude of $(uo) is only 
marginally influenced by the mass. As a result, there is a first order phase transition at any to, and Hc increases with 
TO. The critical curve in the m — fi plane is shown in Figure n, for various Nj.. 
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m 

FIG. 6. The value of the critical chemical potential /ic as a function of the quark mass m, for various A'^ values. The critical 
chemical potential increases with the mass. 

The dependence on the size parameter g is also easy to understand. For small g values, fic approaches zero, since 
the value of the maximum of approaches 1. For large g, u approaches 1 and the value of the maximum grows like 
g^ , so ^.c{g) increases logarithmically. In fact, for zero mass we can compute the limiting curve ^ic{g) at Nt —> oo, in 

a closed form: fic = ^(In A — u) where A=^ + a/i + ^ and u = + /l + This curve and /ic(g) for Nt — 8 
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FIG. 7. The value of the critical chemical potential /ic as a function of the parameter g, for various A'' values and in the 
limiting case. The critical chemical potential increases with g and tends to zero when g vanishes. 

The m and g-dependence of the phase boundaries is nicely illustrated in the next section, by the partition function 
zeros. In particular, for small g the zeros are sitting on a contour around the roots of unity. The contours plotted in 
Figure |ll| are obtained by comparing the value of <i> at u = and some complex /x with that of /i), but with u 
obtained at 3?^ = 0. Even so, the agreement with the line of partition function zeros is very good. 

IV. NUMERICAL RESULTS 
A. Eigenvalue distribution 

It is widely accepted now that the distribution of the eigenvalues in this class of non-hermitian random matrix 
models is connected to the phase structure defined by a partition function where the determinant has been replaced 
by its absolute value While it is likely that the model under consideration is no exception, we will only present 
numerical results here regarding the distribution of the eigenvalues. These are qualitatively similar to lattice eigenvalue 
distributions seen in the literature 1^,^. 

The eigenvalues in the complex m plane are straightforward to define. They are the eigenvalues of the model Dirac 
matrix 

^ V iCt + c^Sat, (g> In, - e^'^S^^ (g) In, 

Because of the chiral structure, they must come in ± pairs. In Figure ^we plot three sets of eigenvalue corresponding 
to ensembles of 10 matrices with Nt — 8 and Ng = 10 for a relatively small value of g — 0.3 and various nonzero fi 
values. We recognize the generic pattern with the eigenvalues originally on the imaginary axis (for /i = 0), which then 
spread out gradually in the complex plane as fj, increases. The qualitative similarity with lattice eigenvalues is stronger 
than for the usual RMM, since here the cloud remains contiguous and only the center is depleted of eigenvalues for 
large ^, similarly to while in the usual RMM the cloud eventually breaks up into two distinct pieces which are 
pushed away along the real axis Q . 
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FIG. 8. Eigenvalues in the mass plane for A^'t — 8, Ns — 10, g = 0.3, /j, = 0.20, 0.40, and 0.80 respectively, from left to right. 
The initial line of eigenvalues along the imaginary axis widens in the real direction. For combinations of large fi and small g 
the eigenvalues are grouped around Nt points located symmetrically inside the ellipsoidal cloud. 

The eigenvalues corresponding to the fj, dependence are perhaps more interesting from a practical point of view. The 
natural variable is the fugacity ^ = e^. We are interested in eigenvalues in the complex fugacity plane, corresponding 
more closely to the ones calculated in the Glasgow method Q in lattice simulations. We seek the complex values of 
the fugacity which cancel the fermion determinant D + m for a given random matrix 'configuration' C. We apply 
Gibbs' trick ||^ to construct a matrix whose eigenvalues are explicitly the zeros in the fugacity plane. Let us define 
the matrices 



m 






m J 























G = 
V = 

so that the fermion matrix is written D + m = G + £,V — ^^IV. It is straightforward to check that 

det{D + m) = det(rV)det(^^+|®^^ ^ ) . (15) 

where we again used the formula for reducing the size of the determinant by 2. Hence we are interested in the 
eigenvalues of the propagator matrix 



^ V J \lNt J \ V ^ 

If the size g of the random matrix is small, the ^-eigenvalues are just the Nt order roots of unity. As g increases, they 
are located in a cloud around those values. The size of the cloud increases with g, to the point where the eigenvalues 
merge into a ring around zero. In figure we plot only the eigenvalues larger than 1 in absolute value (since if ^ is an 
eigenvalue, so is l/<^). 
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FIG. 9. Eigenvalues in ^fe the fugacity plane for Nt = S, Ns = 10, m = and g = 0.15, 0.20, and 0.30 respectively, from left 
to right. Only the ones with < 1 are plotted. For small g the eigenvalues are concentrated around the roots of unity and 
merge into a circular band as g increases. 
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FIG. 10. Eigenvalues in the fugacity plane for Nt = 8, Ns = 10, g = 0.30 and m = 0.0, 0.1, and 0.2 respectively, from left 
to right. Only the ones with l^fcj < 1 are plotted. There is gap around 1 which increases with m. 



B. Partition function zeros 



The zeros of the partition function in exactly solvable models provide an interesting link between finite size and 
continuum properties. In the present case one can easily calculate them by expanding in powers of or m. 

It is a well known phenomenon which we have also discussed in some technical detail in that the zeros trace 
out the analytic continuation of the phase boundaries into the complex plane of the respective model parameter. 

This is illustrated in Figure In], where we show the partition function zeros in the fugacity plane at zero mass, for the 
same Ns = 10 and two different Nt values, 8 and 40. The smaller value corresponds to Glasgow simulations Q, where 
they used eight timeslices. Notice that the zeros are significantly away from a circle, which would be the continuum 
result. In the same figures we plot the analytical line corresponding to the transition in the complex fugacity plane, 
for the same Nt values, but using the saddle point approximation in A^^- 
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FIG. 11. Partition function zeros and the analytical curve that separates the two saddle points in the complex [i plane. 
The curve within the unit circle is just the 'mirror image' of the one outside, since everything is symmetric under /i ^ l//x. In 
the limits of large Nt, the curves become circles = p^c and 1/ ^ic ■ The partition function zeros trace out nicely the saddle 

point result, even for a moderate Ns = 10. 



In Figures |lj and |l^ we plot partition function zeros for Nt — 8 and various values of the quark mass m and the 
size parameter g. As the mass is increased, the two lines of zeros move away from the unit circle. For small values of 
the size parameter g, the zeros sit in loops around the Nt order complex roots of unity. 




FIG. 12. Partition function zeros in the complex fugacity plane, for Nt = S and Ns — 10, for various mass values. 
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FIG. 13. Partition function zeros in the complex fugacity plane, for Nt = 8 and A'^s = 10, for various values of the size 
parameter g. 

V. CONCLUSIONS 

We have constructed a random matrix model for the chiral phase transition at finite chemical potential with a /i 
dependence more similar to that in the lattice implementation of the problem. Our approach may be interpreted as 
ensuring time translation invariance in the large Nt limit. We analyzed the model in the saddle point approximation 
with respect to the number Ng of points per timeslice (or number of instanton modes). Similarly to the usual RMM, 
our model features a first-order phase transition driven by the chemical potential /i. Both the particle number (n) and 
the chiral condensate have a finite jump at this transition. The finite jump, hence the first order character, survive 
when the quark mass is increased. 

We calculated the eigenvalues of the Dirac operator in both the complex quark mass plane and the complex plane 
of the fugacity, e'^. The distributions of the eigenvalues show a strong qualitative similarity to those encountered in 
lattice simulations. We also calculated the zeros of the partition function in the fugacity plane for various parameter 
values and found that their location is well approximated by the phase transition curve obtained from the saddle 
point analysis of the effective partition function. 

The most remarkable property of the model is the ^-independence of the partition function in the low density phase, 
when the number of timeslices Nt is large. As a consequence, the number density is identically zero in this phase. 
This property is different from what is seen in the traditional RMM , where there was a negative number density 
in the low density phase, which coincided with a suppression of the partition function. Our previous conclusion on 
the mechanism that made the Glasgow method exponentially slow points precisely towards this suppression as the 
ultimate source of trouble 0. 

For small Nt values the suppression of the partition function is significant in this model, too. However, it disappears 
quickly as Nt increases, together with the fi dependence in the broken phase, as one would expect in continuum QCD. 
The Glasgow lattice simulations which have been reported Q employed a small number of timeslices, corresponding 
to Nt = 8, Nl = 4 in our model. For this value, we find a significant suppression of the partition function. The 
suppression is further enhanced by the number of spatial points, corresponding to our Ng, which is of order Nf in a 
traditional lattice simulation. This leads to the exciting possibility that lattice (Glasgow) simulations of QCD at finite 
fi might be plagued by the same problem of a spuriously suppressed partition function, as one finds in the RMM. The 
suppression is an artifact of time discretization, and could in principle be cured by using a large number of timeslices 
and a small spatial lattice. 

The main argument against excessive enthusiasm is the fact that in the present model, similarly to the usual RMM, 
we have a major cancellation due to averaging. Averaging over the phases of the determinant achieves a cancellation 
in the broken or low density phase, namely from the magnitude of the typical determinant, of the order of exp(2/i7V), 
a large number, to that of the same determinant at /i = 0, which is of the order 1. This cancellation appears then to 
be similar to the one seen in the usual RMM. However, in that case the determinant at /i 7^ was not significantly 
larger than 1, and the cancellation due to phase averaging led to a small number, exp{—fpN). That cancellation 
coincides with the suppression of the partition function at 7^ with respect to its value at /x = 0. As we have 
learned from the analysis of the Glasgow method applied to that model, a partition function which is suppressed for 
some parameter values is very hard to calculate because of the huge precision required. 

One should make a distinction between the cancellation brought about by phase averaging, i.e., from < \det(p)\ > 
to < det(fi) >, and the suppression of the partition function, i.e., from < det{0) > to < det{fi) >. A partition 
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function like the present one, which is the result of large cancellations during averaging, but does not have significant 
suppression, in not a priori incalculable from simulations, and might be amenable to a calculable form using an 
appropriate resummation [T^ ]. 

Our previous study of Glasgow averaging in the usual RMM indicates the suppression of the partition function, 
Z(/i)/Z(0) << 1 as the source of trouble. The pattern of false zeros was not driven by the competition between 
the true partition function and the Stephanov phase (defined by the absolute value of the determinant), but rather 
between the former and a phase resulting from spoiling the cancellation between various terms of the polynomial 
expansion of the true partition function. It is probably too early at this point to encourage a new set of lattice 
Glasgow simulations. However, a reexamination of the Glasgow method in the context of this model might lead to 
encouraging results. 

Our model is similar to those arising from strong-coupling analysis of lattice QCD at finite fj, but it has the 
simplicity of random matrix models. There is a straightforward connection with RMM models with many Matsubara 
frequencies [|rT| . This can be seen by diagonalizing the hopping matrix in the original form of the partition function. 
Extracting a T — /i phase diagram from this model is also possible. Introducing a dimensionful lattice spacing a leads 
to identifying the quantity aNt with the inverse temperature. However, a more careful analysis of the implementation 
of temperature in this model should precede that. 

In summary, we studied a simple model which shows a clear connection between time translation invariance and 
the physical requirement of having a zero baryon number density in the low density phase of QCD. The careful 
implementation of time dependence seems to have the potential to cure one of the artifacts that make traditional 
lattice simulations at finite density virtually impossible. 

It would be interesting to see how the introduction of an explicit temperature dependence in this model would 
change the /x — T phase diagram that emerges from the usual random matrix model. Also, it would be very interesting 
to see how the exciting results obtained in the RMM for color superconductivity would be influenced by our 
method of introducing the chemical potential. 

This work has been supported in part by a grant from the U.S. National Science Foundation. Thanks are due to 
James Osborn, Dirk Rischke, Dominique Toublan, and Jac Verbaarschot for many useful discussions. Ralph Amado, 
Dirk Rischke and Jac Verbaarschot are thanked for a critical reading of the manuscript. I am grateful to the Nuclear 
Theory Group at BNL for their hospitality during part of the time the paper was written. 



APPENDIX A: DERIVATION OF det{QNiA,X)) 



We wish to study matrices of the form 



QLia,x) 



fax 
V X ... 



X-' \ 



(Al) 



which occur in the dual representation of our partition function. This result can also be found in p4| . We derive it 
here for completeness. The size of the matrix N is an arbitrary positive integer. There are no special restrictions on 
a and x. This determinant is We will calculate the determinant by induction. First, we define an auxiliary matrix, 



Q.Nia,x) 



fa X ... \ 



V 



(A2) 



The respective determinants Qn and Qn verify 



(A3) 
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which is obtained as follows. First expand Qn by its first row, which leads to three sub-determinants, one of them 
Qn-i- Expand again the remaining two by their first column, which leads to two matrices Qn-2 and triangular 
matrices, which give the terms. The same trick applied to Qjv yields a recursion relation, 

Qn = aQjv-i-Qjv-2 . (A4) 



We have Q2 = a? — \ and Q3 = — 2a, which leads to Qo = 1 and Qi = a. The recursion formula (A4) can be solved 
explicitly. First, notice that the recursion is homogeneous so if two distinct series {qk} and {f?^} verify the recursion 
relation, so will any linear combination of the two series. Of course, there may be at most two linearly independent 
solutions, since the first two terms of a given solution completely define the whole solution. We seek the solution in 
the form of a power series, qk = A^'^. The base ^ must then verify the characteristic equation 

f -a^+l = . (A5) 

Consider the two roots of ([A5|), Ai2(a) = ^ (a ± \/a^~^). They are each others' inverse, Ai(a) — j^^, and their 

sum is a. We can then write Qk — aiAj;' + a2Af , where the coefhcient-a ai2 are determined from the first two terms 
in the series, as 

ai = T p ; "2 = T ^ ■ (A6) 

The result Qn = ^ A1-A2 — finally leads to the following expression for the original determinant: 

= A^ + A-^-(-)^(x^+x-^) . (A7) 



where for simplicity we defined A = Ai = i(a± ~ 4) — j- 
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